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By designing a proper unitary operator U, we synthesize NMR 
analogues of Einstein-Podolsky-Rosen states (pseudo-EPR states) using 
generalized Grover's algorithm on a nuclear magnetic resonance (NMR) 
quantum computer. Experiments also demonstrate generalized Grover's 
algorithm for the case in which there are multiple marked states. 
PACS number(s):03.67 

Since the quantum searching algorithm was first proposed by Grover [1], 
several generalizations of the original algorithm have been developed [2]- 
[4]. One of the generalized algorithms can be posed as follows. For a sys- 
tem with a total of N basis states, a composite operator Q is denned as 
Q = — I s U~ 1 ItU. U can be almost any valid quantum mechanical uni- 
tary operator. I s is denned as I s = I — 2\s >< s\, where I denotes unit 
matrix. \s > denotes a predefined basis state that is used as an initial 
state in our experiments. I s is a diagonal matrix with all diagonal elements 
equal to 1 except the ssth elements which are -1. Similarly, I t can be writ- 
ten as It = I — Z^2|£ >< t\, where \t > denote the marked states, and 
there are r marked states. For any \t >, I t \t >— — \t >. uis defined as 
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u = yY,t \Uts\ 2 , where U ts =< t\U\s >. It has been proved that tt/Au ap- 
plications of Q transform the system from \s > into a superposition denoted 
as \ip >= \ Y*tUtsU~ l \t > if u <C 1. Through introducing an ancilla qubit 
and choosing a proper U, Grover proposed a theoretical scheme to synthe- 
size a specified quantum superposition on N states in 0(y/~N) steps using the 
generalized algorithm [5]. Nevertheless, we find that some useful superpo- 
sitions, such as NMR analogues of Einstein-Podolsky-Rosen states (pseudo- 
EPR states) [6] -[8], can be synthesized using the algorithm without an ancilla 
qubit. Such superpositions can be represented as \ip au >= ^X^e^*|£ > (su 
for superposition), where <f> t denote the phases of \t >. By designing a proper 
U, we make \U ts \ identical, and let jj^j = e 1 ^, so that \ip > can be repre- 
sented as \ip >= ^ X^e l< ^[/ _1 |i >. After the application of U, the system 
lies in \ip su >, where an irrelevant overall phase factor can be ignored. 

In our previous work, we have realized generalized Grover's searching 
algorithm for the case in which there is one marked state on a two-qubit 
NMR quantum computer [9]. In this paper, we will synthesize pseudo-EPR 
states using the algorithm. 

Our experiments use a sample of carbon-13 labelled chloroform dissolved 
in d6-acetone. Data are taken at room temperature with a Bruker DRX 500 
MHz spectrometer. The resonance frequencies v\ = 125.76 MHz for 13 C, and 
1/2 = 500.13 MHz for 1 H. The coupling constant J is measured to be 215 
Hz. If the magnetic field is along z-axis, by setting h — 1, the Hamitonian 
of this system is represented as 

H = -2KI/JI - 2-kv 2 II + 2nJllll (1) 

where 1* [k = 1,2) are the matrices for ^-component of the angular momen- 
tum of the spins [10]. In the rotating frame of spin k, the evolution caused 
by a radio-frequency (rf) pulse on resonance along x or — y-axis is denoted as 
X k {ip k ) = e lLpkI * or Y k (-(p k ) = e~ t,fikI v , where y? fe = B x ^ k t p with k specifying 
the affected spin. B±, j k and t p represent the strength of rf pulse, gyromag- 
netic ratio and the width of rf pulse, respectively. The pulse used above is 
denoted as \(p\ k x or [— The coupled-spin evolution is denoted as 



where r is evolution time. The initial pseudo-pure state 



\s >= | T>i I T>5 
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is prepared by using spatial averaging [11], where | \>k denotes the state of 
spin k. For convenience, the notation | |>! | f > 2 is simplified as | tt>- The 
basis states are arrayed as | Tt>> I tl>> I lt>> I ll>- Pseudo-EPR states 
are denoted as |V>i >= (| TT> +| ll>)/v / 2, |^ 2 >= (I TT> -| ll>)/v / 2, 
IV* >= (I T!> +| IT>)/v / 2, and |^> 4 >= (| U> -I IT>)/v / 2. EPR (or 
pseudo-EPR) states are very useful in quantum information and have be 
implemented in experiments [12] [13]. We will synthesize pseudo-EPR states 
using generalized Grover's algorithm. 

U is chosen as U = Y 1 (<f 1 )Y2((p 2 ) represented as 



U = 



( C\C 2 C\S 2 S\C2 S 1 S 2 \ 

-CiS 2 C X C 2 ~S 1 S 2 SiC 2 

-SlC 2 ~S 1 S 2 CiC 2 C 1 S 2 

V SiS 2 SiC 2 -CiS 2 CiC 2 ) 
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where Ck = cos(<^fc/2), Sk = sm(^/2). According to the first column of 
U, we design (p 1 and (p 2 for synthesizing pseudo-EPR states. U is chosen 
as V x = n(f)y 2 (f ), C/ 2 = Yi(f)y 2 (-f ), C/ 3 = y(f)y 2 (f), and u 4 = 
Yi(— I )y 2 (f ) for >, |^2 >, 1^3 >, and |^4 >, respectively. One can prove 
that Uj = \Jj2t \Ujts\ 2 = 1/2, where j=l,2,3 or 4. The following discussion 
will show the condition u <C 1 is not necessary. If the system starts with 
other basis states, U is chosen as other forms. For example, if I s — | Tl>, 
according to the second column of U, U is chosen as U\ = Yi(— f)y(f), 
U 2 = y(f )y 2 (f ), U 3 = y(f )y 2 (-f ), and U 4 = y(f )y 2 (f ). If C/ is chosen 
as fA = ^(^(f ), C/ 2 = X 1 (*)Y 2 (-%), U 3 = X 1 (5)y 2 (f), and C/ 4 = 
A"i(f)y 2 (-f), pseudo-entangled states (| |T> -i| ll>)/v / 2, (| TT> +i\ ii> 
)/V2, (| U> -i| IT>)/v / 2, and (| ||> +i| IT>)/v / 2 can be obtained, 
respectively. 
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If \s >— I tt>, h is chosen as I represented as 



/ -1 \ 
10 
10 

V o 001/ 

Because r = N/2, the conditional sign flip operators for | |"f> and 
and for | ||> and | J,f> can be chosen the same form represented as 

/ 1 \ 
0-100 

0-10 

V 1 / 



(5) 
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Q is represented as Qj = —I s Uj 1 hUj for >. The operator Gj is defined 
as Gf ] = U,Q r l, which means that operation Qj is repeated n times, and 



then Uj is applied. It is easy to prove that G^\ Tt>= — |^i >, G^'l 1 
-\ip 2 >, G { 3 } \ tt>= — 1-03 >, and G!p\ TT>= ~\^4 >■ The required n to 
synthesize a target state displays a period of 3. For example, G^\ tt>= 

i^ 3 >,Gi 7) i tr>=-i^ 3 >. 

The following rf and gradient pulse sequence [a] 2 x — [grad] z — [7r/4]J. — 
1/4J— [n]]: 2 — 1/4J- — n]]: 2 — [— n/A\ l y — [grad] z transforms the system from 



the equilibrium state represented as 



Peg = 

to the initial state represented as 



ji/2 + I 2 /2 + IH 



2 

z 1 



(7) 
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which can be used as the pseudo-pure state | |t> [14]. a = arccos (71/272), 
[grad] z denotes gradient pulse along z- axis, and the symbol 1/4 J means that 
the system evolutes under the Hamitonian H for 1/4 J time when pulses are 
switched off. The pulses are applied from left to right, [tt]^,' 2 denotes a 
nonselective pulse (hard pulse). The evolution caused by the pulse sequence 
1/4 J — [n] 1 ^ 2 — 1/AJ — [— tt]1' 2 is equivalent to the coupled-spin evolution 
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[1/2 J] described in Eq.(2) [15]. [ir]l' 2 pulses are applied in pairs each of 
which take opposite phases in order to reduce the error accumulation causes 
by imperfect calibration of [tt] pulses [16]. 

U u U 2 , U 3 and U 4 are realized by [tt/4]J - [37r/4] 2 , [tt/4]J - [-3tt/4] 2 , 
[7r/4]J - [7r/4] 2 , and [— 7r/4]J - [7r/4] 2 , respectively. J t = [l/J], realized by 
1/2 J - [?r]i' 2 - 1/2 J - [-7T& 2 . According to Ref.[17], /„ is realized by 1/4J- 



7T 



1,2 



- 1/4J - [-tt]^ 2 - [-vr/2]^ 2 - [-vr/2]^ 2 - [tt/2]^ 2 . transforms 
the system from the initial state into the corresponding target state. For 
example, G^p transforms the system from p into p^ represented as 



P 3 = Illl + I l y ll ' (9) 

which is equivalent to j -03 > < -0 3 1 . A readout pulse [7r/2]* transforms p 3 
+ Iyly + III " 2 , which is equivalent to 



to 



P3r = 7 



/ 1 1 1 -1 \ 
111-1 
111-1 

V -l -l -l i , 



(10) 



The information on matrix elements (1,3) and (2,4) in Eq.(10) can be directly 
obtained in the carbon spectrum, and the information on elements (1,2) and 
(3,4) can be directly obtained in the proton spectrum. Similarly, when the 
system lies in \ipi >, \ip 2 >, or >, through the readout pulse [7r/2]*, the 
system lies in the state described as 
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Through observing the matrix elements (1,3), (2,4), (1,2) and (3,4) in Eqs.(lO)- 
(13), one can distinguish pseudo-EPR states. 

In experiments, for each target state, the carbon spectrum and proton 
spectrum are recorded in two experiments. For different target states, carbon 
spectra or proton spectra are recorded in an identical fashion. Because the 
absolute phase of an NMR signal is not meaningful, we must use reference 
signals to adjust carbon spectra and proton spectra so that the phases of the 
signals are meaningful [18]. When the system lies in the initial pseudo-pure 
state described as Eq.(8), the readout pulses [7r/2]J and [tt/2]^ transform it 
into states represented as 



Psrl — -. 



and 



Psr2 — T 
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(14) 



(15) 
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respectively. In the carbon spectrum or proton spectrum, there is only one 
NMR peak corresponding to element (1,3) in p sr i or to element (1,2) in p sr2 . 
Through calibrating the phases of the two signals, the two peaks are adjusted 
into absorption shapes which are shown as Fig. 1(a) for carbon spectrum and 
Fig. 1(b) for proton spectrum. The two signals are used as reference signals of 
which phases are recorded to calibrate the phases of signals in other carbon 
spectra and proton spectra, respectively. One should note that the minus 
elements in Eq.(14) and Eq.(15) correspond to the positive peaks in Fig. 1(a) 
and Fig. 1(b). 
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Experiments start with pseudo-pure state | ft>- Gj ' transforms | |t> 
into >. If no readout pulse is applied, the amplitudes of peaks is so 
small that they can be ignored. By applying the spin-selective readout pulse 
[71-/2]^, we obtain carbon spectra shown in Figs. 2 (a), (b), (c), and (d), and 
proton spectra shown in Figs. 2(e), (f), (g), and (h), corresponding to >, 
\ip2 >, \ip3 >, and \ip4 >, respectively. In Fig. 2(a), for example, the right 
and left peaks correspond to the matrix elements (1,3) and (2,4) in Eq.(ll), 
respectively. Similarly, in Fig. 2(e), the two peaks correspond to the matrix 
elements (1,2) and (3,4) in Eq.(ll). The phases of the signals corroborate 
the synthesis of pseudo-EPR states. 

In conclusion, we synthesize pseudo-EPR states using the generalized 
Grover's algorithm by choosing a proper U. Although the ancilla qubit is 
not used, our experimental scheme shows the essential meaning of Grover's 
original idea. The experiments also demonstrate generalized Grover's algo- 
rithm for the case in which there are N/2 marked states. 

This work was partly supported by the National Nature Science Foun- 
dation of China. We are also grateful to Professor Shouyong Pei of Beijing 
Normal University for his helpful discussions on the principle of quantum 
algorithm. 
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Figure Captions 



1. The carbon spectrum (Fig. 1(a)) obtained through selective readout 
pulse for 13 C [7r/2]* and the proton spectrum (Fig.2(b)) obtained through 
selective readout pulse for 1 H [n/2]y when the two-spin system lies in 
pseudo-pure state | Tt>- The amplitude has arbitrary units. The two 
peaks are adjusted into absorption shapes. The two signals are used as 
reference signals to adjust other spectra. 

2. Carbon spectra (shown by the left column) and proton spectra (shown 
by the right column) obtained through [ir/2] y after pseudo-EPR states 
are synthesized. Figs. 2(a)- (d) and Figs. 2(e)- (h) correspond to states 
(I TT> +1 !l>)/v / 2, (I TT> -I ||>)/v / 2, (| U> +| iT>)A/2, and 
(I U> -I IT>)/v / 2, respectively. 

[Figure 1 about here.] 

[Figure 2 about here.] 
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